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Abstract 

We present predictions for the production cross section of t-quark pair production in association 
with a Z boson at the next-to-leading order (NLO) accuracy using matrix elements obtained 
from the HELAC-Oneloop package. We use the subtraction method for computing the radiative 
corrections as implemented in the POWHEG-Box, which was also used in several other computations 
of similar complexity. 
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In the recent years we have been witnessing the NLO revolution in computing QCD 
jet cross sections. Before the millenium, one-loop amphtudes were computed analytically, 
which was a serious bottleneck for computing QCD radiative corrections to processes of 
final states with high multiplicity. Indeed, the most complex computations involved three- 
jet production with hadrons in the initial state [U E] and electron-positron annihilation into 
four-jets [3j. The emergence of unitarity-based, fully numerical approaches to computing 
one-loop amplitudes P-[TT] has changed the scope of possible computations completely. By 
now there are publicly available computer programs for this purpose [121 US], which paved 
the way to the automation of NLO computations [Ti] . 

In this letter we present predictions for the process pp — t- tiZ at the NLO accuracy in 
QCD. This process is interesting from the point of view of measuring the iiZ couplings at 
the LHC p^J, which has not yet been measured directly. This coupling is predicted by the 
Standard Model (SM). Significant deviation from this predicted value could be possible signal 
of new physics beyond the SM. In order to optimize such measurements precise theoretical 
predictions are needed both for the signal and the background processes. 

Our present computations constitute steps in an ongoing project for generating event 
samples for pp ^ t t+X processes, where X is a hard partonic object in the final state 
[TU| [T7] . These event files are stored according to the Les Houches accord [TH] , and can be 
interfaced to standard shower Monte Carlo programs to produce predictions for distributions 
at the hadron level that are exact up to NLO accuracy upon expansion in the strong coupling. 
With such predictions at hand one can optimize the selection cuts for the signal process for 
improved experimental accuracy of the coupling measurements. 

In our project we use the POWHEG-Box [19J, which produces a general framework also to 
NLO computations based on the FKS subtraction scheme pO]. The box requires the relevant 
matrix elements as external input. For this purpose we use the HELAC-NLO package [211 122], 
from which we obtain the matrix elements in a semi-automatic way. Although these tools 
have already been used successfully for many computations, the complexity of the problem 
requires tedious checking and validation procedures to make sure that the event samples 
can be safely used in experimental analyses. Part of this procedure is a careful check of the 
predictions at the NLO accuracy which is the subject of this letter. 

We performed our calculations using the POWHEG-Box which requires the following ingre- 
dients: 
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• Flavor structures of the Born {gg — t- Ztt, qq — )■ Ztt, qq — > Ztt) and real radiation 
{qq — 7- Zttg, gg — > Ztig, qg — )■ Ztiq, gq — )■ Ztiq, qq — )■ Zttg, qg — )■ Ztiq, gq — )• Ztfg) 
subprocesses (g G {m, (i, c, s, b}). 

• The Born-level phase space was obtained by using the invariant mass of the tt-pair 
and four angles. 

• Crossing invariant squared matrix elements with all incoming momenta for the 
Born and the real-emission processes were built using amplitudes obtained from 
HELAC-Oneloop ^ and HELAC-PHEGAS [21], respectively. The matrix elements in the 
physical channels were obtained by crossing. 

• The color matrices for the color-correlated squared matrix elements were taken from 
HELAC-Dipoles [25]. 

• We used the polarization vectors to project the helicity amplitudes to Lorentz basis 
for writing the spin-correlated squared matrix elements. 

With this input POWHEG-Box can be used to perform the necessary integrations numerically. 
In order to ensure the correctness of the computations, we performed comparisons of the 
matrix elements to those obtained from independent sources and consistency checks. 

We compared the tree-level (Born and real-emission) matrix elements to those obtained 
from MADEVENT [26], while the matix element for the virtual correction, including the first 
three terms in the Laurent expansion in the dimensional regular isat ion parameter e was 
compared to that from GoSam [13j in randomly chosen phase space points. In all cases we 
found agreement up to at least 6 digits. We also computed various distributions at leading 
order with both PowHel (=POWHEG-Box+HELAC-NLO) and MADEVENT and found agreement. 

The consistency between real-emission. Born, color-correlated and spin-correlated matrix 
elements was checked by taking the soft- and coUinear limits of the real-emission squared 
matrix elements in all possible kinematically degenerate channels using randomly chosen 
phase space regions. As a further consistency check, we have also implemented the similar 
process pp — )■ tt7, which has also been computed at the NLO accuracy [2HII2S], and we found 
agreement with the predictions of Ref. [29]. The important difference from the point of view 
of the computations lies in the mass of the vector boson. In this respect the tt7-production 
resembles tt + jet -production presented in Ref. [16], while ttZ-production resembles tlH- 
production, presented in Ref. [T7] . 
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FIG. 1. Transverse momentum distribution of the Z-boson from Ref. [30j (LMMP) and from our 
calculation (PowHel). The lower panel shows our prediction for the NLO iT- factor compared to 
that of LMMP. 



For ttZ hadropro duct ion the fully inclusive cross section and the distribution of the 
transverse momentum of the Z-boson is the only published prediction at NLO accuracy. 
In Ref. [30] the prediction was made with the default scale choice /i(= = \ip) = /io(= 
rrit + mz/2). With this scale and using the same input parameters we find agreement 
for the LO prediction (0.808 pb), including the dependence on the scale fi. However, at 
NLO accuracy our prediction is slightly larger, (1.121±0.002) pb compared to 1.09pb [31] 
in Ref. [SD], leading to an inclusive fCinc = 1-39 (instead of the i^'inc = 1-35 in Ref. [50]). 

For the p^^^-distribution we checked our LO prediction with MADEVENT, and found agree- 
ment as shown in Fig.[l| but not with the LO predictions of Ref. f3ni|. As a result, the 
predictions of the two computations disagree also at NLO. In order to separate the effect of 
the differences at LO, in the lower panel we plot the i^'-factor of our computation against 
the -R'nlo of Ref. [30]. The two predictions agree within the uncertainties of the numerical 
integrations except in two bins. 

The almost constant value of the f^-factor for the p^-distribution of the Z-boson made 
the authors of Ref. [^ expect that the shape of many other kinematic distributions will 
also be approximately unchanged by NLO corrections. Our findings do not support this 
expectation. As an example, we show the p^-distribution of the t-quark in Fig.[2j 

Next we turn our attention to make predictions for ttZ hadroproduction at the current 
LHC The parameters of our calculation were the following: ^/s = 7TeV, CTEQ6.6M PDF 
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FIG. 2. Transverse momentum distribution of the t-quark. The lower panel shows our prediction 
for the NLO JT-factor compared to a constant value i^mc = 1-39. 



set [32] from LHAPDF, 2-loop running a^, with A^^^ ^ 226 MeV, rrit = 172.9 GeV, mz = 
91.1876 GeV, mw = 80.399 GeV, Gp = 1.16639 • 10"^ GeV'^ the renormahzation and 
factorization scales were chosen equal to fiQ = rrit + mzf^. 

In Fig.|3] we show the transverse momentum distributions of the Z-boson and the t- 
quark. We find that the i^'-factor for the p^^^-distribution is somewhat smaller and less 
uniform than that at 14TeV. The bands correspond to the simultaneous variation of the 
renormahzation and factorization scales between one half and two times the default value. 
We observe a significant reduction of the scale dependence as going from LO to NLO, also 
seen for other distributions. The extra radiation, present in the NLO computation, makes 
the spectra softer as expected. This softening decreases the i^'-factor below one for very 
large values of transverse momentum of the t-quark. Around the value of p± where the 
width of the NLO band shrinks to zero accidentally, and the band is unlikely to represent 
the effect of the missing higher order contributions. 

In Fig.|4]we show the rapidity distributions of the Z-boson and the t-quark. The correc- 
tions are clearly not uniform, but increase towards large rapidities, making the predictions 
less reliable in those regions. For central rapidities the corrections are moderate and the 
scale uncertainty decreases significantly. 

Finally, Fig.|5] shows distributions of the separation in the rapidity-azimuthal angle plane 
of the tt-pair and also of the t-quark Z-boson. For the tt-pair the corrections are large 
except for the back-to-back configuration, while for the tZ-separation the corrections are 

5 



2 
1 

CD 

^ 2 
^ 10-' 



2 

10-^ 
1.5 

1.25 

1.0 



1 1 1 1 1 1 1 1 1 

■ ~\- 

=- \/s = 7TeV L 

: CTEQ6.6M 

- m, = 172.9GeV 

■ mz = 91.1876GeV 

;io = mt + m,z/2 


1 1 1 1 1 1 1 1 1 1 

— LO -= 

NLO E 

/j(i/2 < ii< 2/1(1 2 

-— — 

1 1 1 1 1 1 i iT^ 


;i 1 1 1 1 1 1 1 1 




llllll 


= 



1 FT 



50 100 150 200 250 300 350 400 450 500 

P±,z [GeV] 



> 
O 

£ 10- 



5 - 

2 - 

I 



I M ' M M M M 

— LO 

NLO 

'~mtm— //ii/2 < /I < 2/io 



i = 7TeV 
5 : CTEQ6.6H 

mt = 172.9GeV 
mz = 9L1876GeV 
/.lo = m + 171^/2 



2 - 



O 



^5 ^ 



1.25 E- 
1.0 




I I I I I I I I I i : 




50 100 150 200 250 300 350 400 450 500 

P±,t [GeV] 



FIG. 3. Transverse momentum distributions of the Z-boson (left) and the t-quark (right). The 
lower panels show our prediction for the NLO ii'-factors compared to the constant value K^^c = 1-39 
found at 14TeV and default scale fiQ. 
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FIG. 4. Rapidity distributions of the Z-boson (left) and the t-quark (right). The lower panels 
show our prediction for the NLO i^- factors compared to the constant value i^mc = 1-39 found at 
14TeV and default scale //q- 



moderate except for large rapidities. 

In this letter we studied the hadroproduction of a tt-pair in association with a Z boson. 
This process is of interest for measuring the ttZ-coupling directly at the LHC. Our predic- 
tions slightly disagree with the results of the only published computation. The difference 
most likely originates from the lower numerical precision of the integrations in Ref. [30] . 

We produced predictions for the LHC. In general one can conclude that the NLO pre- 
dictions are sizeable, but not alarmingly large, lie in the 30-40 % range except towards the 
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FIG. 5. Distributions of the i?-separation between the tt-pair (left) and between the t-quark and 
Z-boson (right). The lower panels show our prediction for the NLO i^-factors compared to the 
constant value i^mc = 1-39 found at 14TeV and default scale //q- 



edges of the kinematically available regions, where the i^T-factor can grow above 1.5. As a 
result the dependence on the renormalization and factorization scales decreases significantly 
and the theoretical prediction becomes fairly reliable. 

In order that our predictions become useful for the measurement of the ttZ-couplings, 
matching with shower Monte Carlo programs is desirable, which we consider in a separate 
publication. With such a matching we can let the t-quarks decay and study the optimization 
of the experimental cuts realistically. 
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